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iV-Body Nuclear Forces at Short Distances in Holographic QCD 
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We provide a calculation of iV-body (N > 3) nucleon interactions at short distances in holographic 
QCD. In the Sakai-Sugimoto model of large N c massless QCD, N baryons are described by N Yang- 
Mills instantons in 5 spacetime dimensions. We compute a classical short distance interaction 
hamiltonian for N 'tHooft instantons. This corresponds to N baryons sharing identical classical 
spins and isospins. We find that genuine iV-body nuclear forces turn out to vanish for N > 3, at 
the leading order. This suggests that classical iV-body forces are always suppressed compared with 
2-body forces. 



- Introduction. Recent developments in computa- 
tional nuclear physics reveals that the microscopic de- 
scription of nucleus in terms of nucleon degrees of free- 
dom requires three-nucleon interactions. In fact, al- 
though the nuclear three-body interaction is weaker than 
the two-body interaction, the binding energies of light nu- 
clei [l| and the saturation density of nuclear matter 
cannot be understood without taking into account the 
three-body terms. This is due to a large cancellation of 
the kinetic energy and two-body attraction. The main 
component of the three-body interaction is associated 
with two-pion exchange, such as the Fujita-Miyazawa 
force [3J. However, in addition to this, a repulsive three- 
body interaction of short range is required for quanti- 
tative description of nuclear systems [l], The short- 
range three-nucleon interaction, which is assumed to be 
spin-isospin independent in many cases, is important for 
determination of the nuclear equation of state at high 
density @,[|. 

The two-body interactions adopted in those many- 
body calculations are determined by the phase-shift anal- 
ysis of nucleon-nucleon scattering data. However, much 
less information is available for the iV-body forces (N > 
3). Of course, we know that, in principle, the nuclear 
properties should be derived from QCD [6(. However, 
QCD is strongly coupled at the nuclear energy scale, 
which leads to a huge gap between QCD and nuclear 
many-body problems. 

A recent progress in string theory can bridge this gap, 
analytically. It is called holographic QCD, an application 
of of gauge/string duality Q to strongly coupled QCD. 
We apply the holographic QCD to iV-body nuclear force 
(JV>3). 

In holographic QCD, one of the most successful D- 
brane models is Sakai-Sugimoto model (SS model) fl, [§]. 
The theory, which is a U(Nf) Yang-Mills-Chern-Simons 
(YM-CS) theory in a warped 5-dimensional space-time, 
was conjectured to be dual to low energy massless QCD 
with Nf flavors, in the large iV c and large A limits (A = 
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N c gQ CD is a 'tHooft coupling of QCD). Modes of the 
gauge fields correspond to meson degrees of freedom and 
this model reproduces surprisingly well various expected 
features of hadrons, incorporating very nicely the nature 
of chiral Lagrangians. 

Baryons are identified with soliton solutions localized 
in the spatial 4-dimensions Q. This is quite analogous 
to that in pion effective theory, baryons are identified 
with Skyrmions. Quantization of a single soliton in the 
SS model [TH EU gives baryon spectra, and also chiral 
properties such as charge radii and magnetic moments 
[T3 | (for other approaches to baryons, see [Hj]). Meson- 
baryon-baryon couplings [12J give a basis of a 2-body 
nuclear force at long distances, a la one-meson-exchange 
picture. Short-distance nucleon-nucleon forces were com- 
puted [14( | , which generates a repulsive core with analytic 
formula for potentials in the large N c limit. A key is 
that the warping can be absorbed into the rescaling of 
the YM-CS theory and brings the string scale to QCD 
scale. Furthermore, when two solitons are close to each 
other, the warping factor is almost constant, therefore 
the effects of the curved geometry can be ignored so that 
an exact two-soliton solution is available. 

In this letter, we compute iV-body nuclear forces for 
arbitrary N, with exact iV-instanton solutions, general- 
izing the method in Ref. [l4[ . The exact treatment is in 
contrast to the Skyrmion and other chiral soliton mod- 
els, in which multi-soliton solutions are quite difficult to 
obtain. 

— Nuclear Force at Short Range. Baryons, includ- 
ing nucleons, are identified with solitonic solutions in the 
SS model @ , and we provide a brief review of the con- 
struction of the soliton and the 2-body nuclear force at 
short range for two flavors (Nf = 2) computed in [Til ]. 

The following rescaling of the coordinates can al- 
low one to understand the system as a 1/A perturbation 
around a flat space, which is suitable for studying the 
instanton solution: 

~m = x i/2 x m ( M = 1,2,3,4), x° = x°, (1) 

and accordingly Ao(t, x) = Ao(t, x) and Aj\i(t, x) — 
\-i/ 2 A]M;(t,x). In the following we omit the tilde for 
simplicity. In these new variables, there are essentially 
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two deviations from the YM theory in the flat space, at 
the leading order in 1/A expansion: (i) the effect of the 
CS term, and (ii) the effect of the space weakly curved 
along the x A direction. The additional hamiltonians are 
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respectively, with a = l/(2167r 3 ). Note that we work 
in the unit Mkk = 1, where Mkk is the unique scale 
parameter appearing in the model, and it can be fixed 
by fitting the p meson mass, giving Mkk = 949 [MeV]. 

In particular, for a single baryon, the leading order so- 
lution is a single instanton in the flat 4-dimensional space, 
that is BPST instanton [15| . The instanton has moduli 
parameters: the instanton location X , the size p and 
the orientation in SU(2). These hamiltonians induce po- 
tentials in the moduli space of the instanton, and p and 

X 4 prefer particular values classically, (p c \) 2 — §^3^ y§ , 

X^ = 0. For multi-instantons, there appears a poten- 
tial for the moduli representing the distance between the 
instantons, which is in fact the nuclear force. 

In Ref. [lij], this 2-body nuclear force was evaluated 
explicitly. Multi-instanton solutions are available in flat 
space, while in this particular curved space it is diffi- 
cult to find them. However, when instantons are close 
enough to each other, the effect of the curved space 
can be neglected, and as a leading order solution we 
can use the multi-instanton solutions in the flat space. 
Therefore, the distance ry between the i-th and the j-th 
nucleons allowed in this approximation is |ry| < M^ K 
{\rij\ < A 1 / 2 M^) in the original (rescaled) coordinates. 
Thus we probe only the short range for the nuclear force. 

The construction of the two-instanton solution owes to 
the renowned ADHM (Atiyah-Drinfeld-Hitchin-Manin) 
method 0, [nl- The moduli parameters of generic N 
instanton solutions are completely encoded in the real 
N x N matrix function L(x; X, ■ ■ ■). The Osborn's for- 
mula fl8j tells us the instanton density 



tr(F M jv) 2 = Q 2 logdetL 
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where □ = <9m<9m- Using this expression, the equation 
of motion for the U(l) part of the gauge field which is 
sourced by the instanton density is solved as [l(| 
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With this explicit dependence on the instanton moduli 
parameters in L, one can compute the hamiltonians ([2]) 
and ([3]) as functions of them. Then, the expectation value 
of the hamiltonians for given baryon states (the wave 
functions are written by the moduli parameters) gives 
the nuclear force at short range ■ 

There is the third contribution to the additional hamil- 
tonians, if kjn, which is present only in multi-instanton 



case. This comes from the metric of the instanton mod- 
uli space. In Ref. [l4| , it was shown that it is higher order 
in l/N c compared to the other two hamiltonians ((2]) and 
©, so we need not compute it in this paper. 

- 3-body Nuclear Force. The 2-body nuclear force 
computed in Ref. [14j is for generic spin/isospin compo- 
nents. But since explicit generic N instanton solution is 
not available, we consider a special solution called 'tHooft 
instanton which has 5N — 3 moduli parameters (while 
generic instanton solution has 8N — 3 moduli parame- 
ters). It is important to notice that once we restrict our 
moduli space by hand like this, we cannot get the generic 
expression for the nuclear force for given baryon states. 
Instead, what we will obtain is a classical analogue of the 
nuclear force. 

The moduli parameters of the 'tHooft instantons are 
only the size pi and the location Xf* of each instanton 
(i = 1, 2, ■ • ■ , N). The missing parameters, the orienta- 
tions of the instantons in SU(2), are responsible for the 
spin/isospin wave functions of the baryons. Thus our 
analysis with the 'tHooft instantons is restricted to "clas- 
sical" baryons, where all the spin/isospins of the baryons 
are identical classically. 

First, let us show that H^ t ^ given in Eq. ^ is irrel- 
evant to the three-body nuclear forces. We can use the 
generic formula obtained in Appendix C of Ref. 1J| , 



A' 



d 4 x (x 4 ) 2 tr(FMAr) 2 = Stt 2 V (2(X 4 ) 2 + pf) , (6) 



for the N 'tHooft instantons. The expression consists 
of just a sum of each instanton sector, which means that 
there is no term involving the inter-nucleon distance, that 
is, no contribution to the nuclear force. Therefore, we 
compute the other hamiltonian @ in this paper. (The 
contribution from -ffkin is suppressed as in the case of the 
2 instantons (l9j.) 

In this section, we concentrate on the case for N = 3, 
i.e. the 3-body force. 

For three 'tHooft instantons, which correspond to nu- 
cleons sharing classically identical spins/isospins, we have 
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where we omit the index M and denote (x M ) 2 by x 2 . 
Then the Osborn's formula becomes particularly simple, 

logdetL = ^log( a; -A l ) 2 +log/, (8) 

i 

2 

with / = 1 + J2i ( x _jf-)2 ■ This gives the U(l) gauge field 
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The first term in Eq. © is a self-energy which was 
already computed, and / is a harmonic function, i.e., 
□/ = 0. Thus, all we need to evaluate is only the last 
term in Eq. ©, {d M f) 2 /f. 

For three instantons, we can expand the expression for 
(x - X x ) 2 < (x - X 2 ) 2 ,(x 
can approximate (x — X2) 2 ~ (X\ — X2) 2 
a similar expression for X13. Furthermore, for simplicity 
we put Xi 1 — 0. Then, the expansion is 
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We like to compute the potential $2$. As seen from 
the expression for A , the leading term of the nontrivial 



three-body force has the form 
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This means that, 



we should expect the nuclear force appearing in proper 
to the 3-body is, at the leading order, 
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is fixed to be (p) 4 by a dimensional analysis. The expec- 
tation value of this (p) 4 at the leading order in large N c 
is given by the classical value given before. Then, rescal- 
ing the coordinates back as X12 — * \ x ^ 2 Xi2 and write it 
as the 3-dimensional inter-nucleon distance 7*12 since we 
substitute the classical value Xf = 0, we obtain, at the 
leading order in l/N c , 
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again • • • represents term obtained by permutation for 
the indices 1, 2, 3. 

Note that we are working in a regime A -1 / 2 <c 
-X12.13 *C 1 in the unit A/kk = 1- The natural scale for 
the 2-body force [l4| is 0(N C / XX 2 2 ). So, if we consider a 
natural separation of the nucleons as Xij ~ 1/Mkk> the 
3-body force is suppressed compared to the 2-body force. 
We conclude that the 3-body force at short range is small, 
for baryons carrying classical and equal spin/isospins. 

N-body Nuclear Force. We can easily extend the 
analysis in the previous section to N 'tHooft instantons. 
The result for the leading term vanishes again, as we 
explain briefly below. 

The quantity necessary for computing Aq is 



in the original coordinates. Let us consider only terms of 
this leading form (fTTj) . We obtain 
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Performing the derivatives, and using the following inte- 
gration formulas 
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The expansion analogous to 3-body case implies that the 
leading order of the short-range nuclear force in proper 
to iV-body would be 
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Again, we have rescaled back the coordinates to the orig- 
inal coordinates and • ■ • represents permutation terms. 
In the previous section, we showed that for N = 3 this 
leading contribution vanishes, for the 'tHooft instantons. 
In this section, we prove that for any N this leading con- 
tribution vanishes. 

First, in the integral it is straightforward to obtain 



we find that the right hand side of Eq. (|12j) vanishes. 
Therefore, the leading term of the order 1/(X 2 2 A 2 3 ) van- 
ishes. This means that the expansion starts from the 
next-to-leading order, 
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where ■ ■ • represents terms obtained by permutation for 
the indices 1, 2, 3. Here the dependence on pi (i = 1, 2, 3) 
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Using the formula (|T3|) . this can be easily integrated with 
J d 4 x to give 0. 
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In the hamiltonian (fT7|) . there are additional terms 

coming from the first term in Eq. ([9]) . In the same manner 
this is shown to vanish at the order (fTT|) . Therefore, 
we conclude that the leading order TV-body nuclear force 
(jTTJ) vanishes, for arbitrary N. Note that the 2-body 
nuclear force does not vanish at the leading order, as 
N = 2 computation is exceptional. 

- Summary and Discussions. Using the SS model 
of holographic QCD, we have found that the ,/V-body nu- 
clear force at short range (N > 3) is order of N c (Ar 2 )~ N , 
for nucleons sharing identical classical spin/isospins. 
This is small compared to the 2-body force which is 
0(N c (Xr 2 )" 1 ) in contrast, and it leads to a hierarchy 
of the (N + l)-body / TV-body ratio V iN+1 ^> /V^ ~ 
l/(Ar 2 ) < 1 for N > 3, in the unit M KK = 1. This 
suppression is consistent with our empirical knowledge. 

Effects of the short-range many-body interaction be- 
comes more prominent for higher-density nuclear matter. 
Therefore, for physics of neutron stars and supernovae, 
for instance, properties of TV-body interactions such as 
what is revealed in this letter are important, even if qual- 
itative. 

Our computation is not fully satisfactory since the 
quantum spin/isospin states of each baryon have not 
been incorporated. The wave function of the classical 



spin/isospin is a delta-function of the SU(2) orientational 
moduli of the instantons, thus it is difficult to relate it 
with the quantum spin/isospins. Nevertheless, it is quite 
remarkable that the generic TV-body nuclear force can be 
obtained by analytic computations. The successful per- 
formance of this computation owes in particular to the 
simplicity of the SS model, in contrast to other chiral 
soliton models. 

Furthermore, the theory on which our computations 
of the nuclear force is based is not a phenomenological 
model but a theory which has been "derived" from large 
N c QCD at strong coupling, through the gauge/string 
duality. Therefore, in principle, we can try to address 
theoretically what is different from QCD and what is 
inherited from it. The present computations go beyond 
the limitations of arguments using universality of chiral 
symmetry breaking. 

Although the model at hand is for large N c QCD, these 
two properties of baryons in holographic QCD would be 
sufficiently strong motivations for studying holographic 
QCD and its relation to nuclear physics further. 
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